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We investigate the consequences of space-momentum correlations in quark phase-space distri- 
butions for coalescence processes at the hadronization transition. Thus far it has been proved 
difficult to reconcile such correlations with the empirically observed constituent quark number scal- 
ing (CQNS) at the Relativistic Heavy-Ion Collider (RHIC). To address this problem we combine 
our earlier developed quark recombination model with quark phase-space distributions computed 
from relativistic Langevin simulations in an expanding Quark-Gluon Plasma (QGP). Hadronization 
is based on resonance formation within a Boltzmann equation which recovers thermal equilibrium 
and obeys energy conservation in the quark-coalescence process, while the fireball background is 
adjusted to hydrodynamic simulations of semi-central Au-Au collisions at RHIC. To facilitate the 
applicability of the Langevin process, we focus on strange and charm quarks. Their interactions in 
the QGP are modeled using leading-order perturbative QCD augmented by effective Lagrangians 
£^ , with resonances which smoothly merge into hadronic states formed at T c . The interaction strength 

S ■ is adjusted to reproduce the empirical saturation value for the quark-elliptic flow, ~ 7-8%. The 

resulting cj> and J/xj) elliptic flow recover CQNS over a large range in transverse momentum (pr) 
Q\ j within a few percent. As a function of transverse kinetic energy, both the quark spectra from the 

Langevin simulations and the meson spectra generated via resonance recombination recover CQNS 
from zero to at least 3 GeV. 
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I. INTRODUCTION 



| The mechanism of hadronization, i.e., the conversion of quarks and gluons produced in hadronic or electromagnetic 
. reactions into colorless hadrons, is a non-perturbative problem that is presently not calculable within the theory of 
' strong interactions, Quantum Chromodynamics (QCD). For partons produced at large transverse momentum, p t , 
, the factorization theorem of QCD allows to treat the hadronization process via a so-called fragmentation function, 
which is universal and can, in principle, be determined empirically. At low momentum this scheme is no longer 
applicable and other hadronization mechanisms become relevant. The quark coalescence model (QCM) has provided 
\Q • a phcnomenologically successful framework to understand several non-perturbative features of hadron production in 
hadronic collisions. In elementary (p-N and n-N) reactions, flavor asymmetries in kaon and charmed hadron spectra 
have been associated with the recombination of produced strange and charm quarks with valence (and sea) quarks in 
target and projectile [3, H, In heavy-ion reactions at the Super Proton Synchrotron (SPS) [3] and the Relativistic 
Heavy-Ion Collider (RHIC), recombination of quarks from a thermalized Quark-Gluon Plasma (QGP) @, i, 0, & ©, 



rieavy-ion collider 

[ToL [ill . Il2l . [13 , [p4| gives a simple and intuitive explanation of several unexpected features in the observed hadron 
spectra, most notably the large baryon-to-meson ratio and the rather universal constituent quark number scaling 



■ (CQNS) of the elliptic flow coefficient, V2.h{pr) = n q V2, q (PT /n q ) (n q : number of valence quarks in hadron h, pt'- 
transverse momentum of h). This scaling relation implies the momenta of the coalescing quarks to be collincar, and 
its implementation in QCMs usually restricts their applicability to sufficiently large momenta so that the associated 
non-conservation of energy in the hadron formation process is small. Since at high pt parton fragmentation is 
expected to take over, the typical range of applicability of QCMs is at intermediate momenta, 2 GeV < pr < 6 GeV. 
In our recent work we have suggested a reinterpretation of quark coalescence in terms of hadronic resonance 
formation, implemented via q + q — > M scattering into a Boltzmann equation (AT: meson). Energy conservation 
is obeyed by utilizing hadronic reaction rates, along with detailed balance, based on pertinent spectral functions. 
In addition, we have shown that this approach correctly recovers the thermal equilibrium limit, which enabled a 
more controlled extension of the coalescence mechanism to low pt and to make contact with the phcnomenologically 
successful hydrodynamic description of bulk matter at RHIC. 

Another aspect that has evaded a satisfactory explanation in QCMs at RHIC is the question of space-momentum 
correlations in the underlying (thermal) quark distribution functions (see, e.g., Ref. fi"H for a recent critical review). 
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In hydrodynamic models the elliptic flow of produced hadrons is a collective effect that implies a definite correlation 
between the particle's momentum and its spatial position in the fireball, i.e., a (locally thermalized) fluid cell moving 
into a specific direction preferentially emits hadrons in that same direction. Such a correlation is neglected within 
the so-called "factorized" implementation of the parton V2, q {pr) which does not carry any spatial dependence (and is 
therefore identical regardless of the quark's position inside the fireball). While this approximation straightforwardly 
recovers the empirical constituent-quark number scaling (CQNS) of the hadron elliptic flow, it is at variance with the 
hydrodynamic description of V2 as a collective expansion effect. 

Previous attempts to incorporate space-momentum correlations into QCMs have found the empirically observed 
CQNS to be rather fragile [H, H3, EH • Part of the problem is the construction of a realistic transition from the thermal 
to the kinetic regime of the underlying parton phase-space distribution functions, as characterized by the "saturation" 
(leveling off) of the empirical parton i> 2 at about p t ~ 1 GeV. In Ref. [l6| several elliptic "deformations" of a thermal 
blast- wave parameterization have been considered, motivated by different plausible realizations of v%. While some 
features could be ruled out being incompatible with the empirical CQNS, other assumptions did not spoil the latter. 
In Ref. [HI a reduction of the boost velocity at higher momenta was introduced, entailing a violation of CQNS 
at the 20% level. It therefore seems that purely phenomenological prescriptions of space-momentum correlations 
and associated «2 did not arrive at a conclusive interpretation of the key features underlying the quark distribution 
functions. In Ref. [l7| . based on numerical transport simulations, it was even argued that rather delicate cancellations 
must be at work to obtain CQNS for the thermal components, thus raising doubts on the robustness of the coalescence 
approach. In view of the broad empirical applicability of CQNS across different ccntralities, system sizes and collision 
energies [ljl[2(|, such an interpretation would be difficult to reconcile with experiment. 

In the present paper we adopt a microscopic approach to compute quark distributions in four-dimensional phase 
space (transverse position and momentum) by employing rclativistic Langcvin simulations for strange and charm 
quarks within an expanding thermal QGP background. In a strict sense, the underlying Fokkcr-Planck equation is 
applicable for a diffusive treatment of heavy and/or high-momentum quarks, i.e., in a regime where the momentum 
transfers from the heat bath are small. Our simulations for low-momentum (pt < 1 GeV) strange quarks are thus at 
the boundary of applicability of a Fokker-Planck treatment and may be considered as extrapolations thereof. The 
Langevin approach has the attractive feature that it naturally encodes the transition from a thermal to a kinetic 
regime. In particular, when simulating heavy-quark (HQ) motion in an expanding QGP fireball for non-central 
collisions, this transition reflects itself in a saturation of the elliptic flow [HLH^, a key ingredient to CQNS in light - 
hadron spectra observed at RHIC. It turns out that Langevin simulations preserve the i>2 saturation feature when 
applied to strange quarks (with thermal masses of ^0.5 GeV). We will therefore investigate whether the resulting 
quark distribution functions, evolved to the hadronization transition and injected into our resonance recombination 
approach, allow for a better (microscopic) understanding of space-momentum correlations in the coalescence process. 
In view of the rather delicate dependence of CQNS on these correlations (as discussed above) , a realistic treatment of 
the kinematics in the hadron formation process is mandatory, including energy-momentum conservation, non-collincar 
kinematics, and a well defined equilibrium limit. The recombination approach developed in Ref. fl2| satisfies these 
requirements. In addition to this, the QGP evolution and subsequent hadronization are linked via the ansatz that 
resonances play an essential role in hot QCD matter around T c . This scenario is consistent with effective potential 
models where a non-perturbative description of the strongly coupled QGP (sQGP) is realized via bound [23| and/or 
resonance [13, [H| states of deconfined partons. Recent lattice QCD computations support the picture of various (light 
and strange) hadronic states surviving up to temperatures of ~1.5-2 T c [2a. 1271]. 

Our article is organized as follows. In Sec. [IT] we review our earlier developed model [ljj for resonance hadronization 
based on the Boltzmann equation. In Sec. Mil we elaborate the computation of the phase-space distributions of quarks 
obtained from Langevin simulations of an expanding QGP fireball at RHIC. In Sec. IIVI we discuss the numerical 
results for the i>2 coefficients of <f> and J/ip mesons within our model, and discuss their properties in terms of CQNS 
in both transverse momentum and transverse kinetic energy, Kt- Sec. [V] contains our conclusions. 



II. RECOMBINATION FROM THE BOLTZMANN EQUATION 



Following Ref. 12j, our description of hadronization at the critical temperature, T c , is based on the Boltzmann 
equation using resonance quark-antiquark cross sections to compute meson spectra in terms of underlying anti- /quark 
phase-space distributions, f q .q (baryons could be treated in a similar way, e.g., in a two-step process using subsequent 
quark-quark and quark-diquark interactions; in the present paper, we will focus on mesons). The meson phase-space 
distribution, F^j, is determined by the equation 

d A r 

— +v-V)F M (t,x,p) = F M (t,x,p) + f3(x,p) , (1) 
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where p and x denote three-momentum and position of the meson, M, and v — p/Em(p) {n%, Em{j>)=\J m 2 +P 2 - 
meson mass and energy). The total meson width, T, is assumed to be saturated by the coupling to quark-antiquark 
states, M q + q and taken to be constant, with the factor 7 p = Euip)/ 171 accounting for Lorentz time dilation, see 
also Ref. fllj | . Integrating over the fireball volume leads to the momentum-distribution function of the meson, Jm, 
and the pertinent transport equation 

f M (t,p)= [ d 3 xF M (t,x,p), ^- f M ( t ,p)=-—f M (t,p)+g(p). (2) 
J ot 7 P 

The drift term vanishes upon integration over x since it is a total divergence: v ■ V/ M (i, %,P) = V • [vfAi(t, x,p)}. The 
relation of the gain term, g(p), to the underlying microscopic interaction is given by 

f f d 3 pid 3 p2 f 

9(Pl = J d 3 x(3(x,p) = J 1 2 J d 3 x f q (x,&i) fq{x,p 2 ) <t{s) v ro \{pi,p 2 ) & (3) {p-pi - fii) (3) 

with <t(s) the cross section for the process q + q — > M at ccnter-of-mass (CM) energy squared, s = (p^ + p^) 2 , 
where pf\ are the four-momenta of quark and antiquark. The quark phase-space distribution functions are normalized 
as Nq t q = J d ^2 7 ^yd P fq,q{x,p). Throughout this paper, quarks will be assumed to be zero-width quasi-particles with 
an effective mass m q (which contains both thermal and bare contributions). The classical nature of the Boltzmann 
equation warrants the use of classical distribution functions for all the particles, and we assume zero chemical potentials 
for all quark species. The cross section is approximated by a relativistic Breit-Wigner form, 

a(s) = g 477 (Fm)2 (4) 

k 2 (s — to 2 ) 2 + (r?n) 2 ' 

where g a = gui /(g q gq) is a statistical weight given in terms of the spin (-color) degeneracy, <?m (gq.q), of the meson 
(anti-/quark) , and k denotes the quark three-momentum in the CM frame. With M ±=> q + q being the only channel, 
it follows that Tin = r out = T. Detailed balance requires the same T in the loss term on the right-hand side of Eq. (JTJ) , 
thus ensuring the correct equilibrium limit with r — 1/T the pertinent relaxation time. This formulation conserves 
four-momentum and applies to all resonances M with masses above the qq threshold, i.e., for a positive Q value, 

Q = m - (m q + mq) > 0. (5) 

If the 2 — > 1 channel proceeds too far off-shell, i.e., Q < and V < \Q\ (e.g., for pions), other processes need to 
be considered, e.g., q + q — > M + g (which, in principle, is possible in the present framework by implementing the 
respective cross sections). We note that the majority of the observed pions are believed to emanate from resonance 
decays (p, A, a\ etc.); in addition, hydrodynamic calculations suggest that the elliptic flow in heavy-ion collisions 
at RHIC does not change much after hadronization [28|. We also note that in the absence of a confining interaction 
individual quarks and antiquarks remain a part of the heat bath. 

The equilibrium limit is readily recovered within our approach by imposing the stationarity condition, 

^-f M (t,p)=0. (6) 



Then ^ is immediately solved by 



Q(P) = ^9(P) , (7) 



which represents the large time limit of the Boltzmann equation and is the expression that comes closest to the con- 
ventional QCM approximation. For hadronization times less or comparable to the relaxation time, r, the equilibrium 
limit will not be reached and an explicitly time-dependent solution is in order. We have verified numerically that 
Eq. ([7]) accurately recovers the standard thermal Boltzmann distribution for a meson M at temperature T , if the 
constraint of a positive Q value is satisfied (for negative Q the 2 — > 1 channel is inoperative), as illustrated in Fig. [1] 
for the case of the (j) meson: when using thermal quark input distributions with radial flow (blast- wave model), the 
computed equilibrium expression ||7J) is in excellent agreement with the same blast-wave expression (identical temper- 
ature and flow profile) directly applied at the meson level. This reiterates the close connection between equilibration 
and energy conservation in the approach of Ref. [l2j|, providing a significant improvement of previous QCMs. 
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FIG. 1: (Color online) pr spectra for <j) mesons using the Boltzmann recombination equation in the equilibrium limit, Eq. (0, 
based on blast- wave input distributions for strange quarks (solid line), compared to <f> spectra directly obtained from a blast- 
wave expression with the same fireball parameters (temperature T = 180 MeV, radial expansion surface velocity /3o = 0.55c). 
The 0-meson resonance parameters are — 1.02 GeV, = 50 MeV, and the s-quark mass is m s — 0.45 GeV. 

III. PARTONIC SPECTRUM 

While the blast- wave model provides a convenient description of the thermal component of empirical hadron (and/or 
parton) spectra, the transition to the kinetic, and eventually hard-scattering, regime is more involved, especially with 
regard to phase-space correlations within QCMs, as discussed in the Introduction. In an attempt to generate realistic 
quark input distributions for meson formation processes in the vicinity of T c , we here adopt a Fokker-Planck approach 
for test particles evolving in a thermally expanding QGP background. The latter is parameterized with guidance 
from hydrodynamic models for central and semicentral Au-Au collisions at RHIC, implementing empirical values 
(i.e., adjusted to experiment) for bulk matter properties such as total entropy, radial and elliptic flow. The elliptic 
flow is parameterized in terms of a flow profile, whose direction at each transverse position is perpendicular to confocal 
elliptic isobars. Its magnitude is chosen to increase linearly in distance from the center with an average boundary 
value of /?o = 0.55c at the end of the mixed phase of the fireball evolution, i.e., at the hadronization time. The 
acceleration is adjusted so that the fireball at hadronization is approximately circular, with bulk vi^ q ~ 5.5%. This 
model has been employed before in the context of heavy-quark (HQ), i.e., charm and bottom spectra at RHIC [2l| . 
and good agreement with hydrodynamic simulations has been found p2j for the same HQ diffusion coefficient. HQ 
observables, which at present are semileptonic single-electron decay spectra p9l. [30j. exhibit an unexpectedly large 
suppression and elliptic flow which cannot be understood within perturbativc QCD (pQCD) including both radiative 
and elastic scattering. The key microscopic ingredient in Ref. [2l| are resonant heavy-light quark interactions mediated 
via effective (broad) D and B mesons in the QGP [3l|, inspired by the findings of thermal lattice QCD. In connection 
with a "conventional" coalescence afterburner [3^ | at T c , the predictions for single-electron suppression and vi turned 
out to be in fair agreement with data [29|, [3(| ■ In more recent work [25| , the effective interactions have been replaced 
by in-medium T-matrices based on finite-temperature HQ potentials extracted from lattice QCD; these calculations 
not only confirmed the interaction strength generated by heavy-light resonances in an essentially parameter-free way, 
but also identified prc-hadronic meson and diquark channels as the most relevant ones. This, in turn, provides a 
direct link between two main discoveries at RHIC, namely the strongly interacting nature of the sQGP and quark 
coalescence from a collective partonic source. 

In the present paper, we build upon the above findings by extending the Fokker-Planck approach to strange (s) 
quarks. While its applicability criterion, m t 3> q ~ T (m t = m? + pf : transverse mass, q: momentum transfer in a 
typical scattering), seems to be only marginally satisfied for momenta pt < 1 GeV (at least in the early phases of the 
QGP evolution), we note that most of the fireball evolution occurs for temperatures close to T c . At higher p t (or nit) 
one enters the kinetic regime where the Fokker-Planck treatment becomes reliable again. With these limits properly 
satisfied, one may hope that the Langevin simulations also accomplish a reasonable description of the low-pt regime 
(Pt i$ 1 GeV) for strange quarks, including realistic space-momentum correlations, which is one of the main objectives 
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in our work. It remains to specify the interaction strength of the strange- and charm-quark species. Here we take 
further guidance from phenomenology by requiring that the final quark V2(pt) exhibits the characteristic saturation 
(or maximum) value of ~7.5%. Our baseline interaction for the stochastic Langevin force is clastic pQCD scattering 
with a rather large value of a s = 0.4 (which can be thought of as containing radiative and/or parts of non-perturbative 
contributions). However, additional non-perturbative interactions are necessary to achieve a sufficiently large vi. As 
in Refs. [U, [3l| we associate these with mesonic resonance states with an interaction strength controlled by the 
resonance width (with a larger width implying stronger coupling). For s quarks (m s = 0.45 GeV) the "heavy-light" 
resonances re quir e a width of T S g ~ 0.3 GeV, and Y C q ~ 0.6 GeV for c quarks (m c = 1.5 GeV), which is compatible 
with Refs. [2ll. 1 3 1 j] 1 - This hierarchy is qualitatively consistent with the general expectation that resonance/bound-state 
formation is suppressed with decreasing constituent mass (and also borne out of the microscopic T- matrix calculations 
for c and b quarks in Ref. [25|). Finally, we have to specify the initial quark distributions. For c quarks we use the 
initial spectra as constructed in Ref. [2l| , as to reproduce D-meson and semileptonic electron spectra in p-p and d-Au 
collisions. A similar procedure is adopted for strange quarks: we parameterize the quark spectra as a superposition 
of exponential and power-law spectra in a way that experimental kaon spectra in 200 GeV p-p collisions are properly 
reproduced (using 5-function fragmentation into kaons at half the parent-quark momentum; as usual in QCMs, the role 
of gluons is suppressed). In AA collisions, the exponential "soft" part is then scaled with the number of participants, 
ATpart, aud the power-law "hard" component with the number of collisions, N co \\, for a given ccntrality. With the 
interaction strengths and initial conditions fixed, the quark phase-space distribution in semi-/central Au-Au collisions 
are predicted from the Langevin simulations at the end of the QGP (mixed) phase without further adjustments, and 
serve as an input for the meson formation processes as described in the previous Section. The framework developed 
here, i.e., QGP evolution with resonance rescattering and recombination at T c , will be referred to as a "Resonance 
Recombination Model" (RRM). 

The quark phase-space distributions resulting from the Langevin approach embody strong correlations between 
spatial and momentum variables. E.g., quarks at high p t tend to be located in the outer layers of the fireball with 
a preferential alignment of the momentum and position vector directions (i.e., the quark momentum tends to point 
"outward"). Likewise, the collective (radial and elliptic) flow, which implies a well-defined (hydro- like) correlation 
between the position of the fluid cell and its radial motion, imprints this correlation on the (partially) thermalized 
components of the Langevin-generated quark spectra. We recall again that the often used factorized implementation 
of the «2 coefficient in coalescence models completely ignores these rather elementary dependencies. 

The proper implementation of the differential phase-space information, carried by the quark distribution func- 
tions (which is essential for a realistic discussion of hadronic ^-scaling properties), into the hadronization formalism 
requires a few technical remarks. Since thermalization in the longitudinal direction is somewhat controversial, we 
assume the quark distributions to be homogeneous in the spatial z-coordinate and flat in rapidity. This leaves four 
independent transverse variables for each particle, which we choose in azimuthal form, (p t , (\> v ,r t , 4> r ), corresponding 
to the distribution dN q /d 2 pt d 2 r t . This 4-D phase-space is then divided into finite bins (with a maximum value of 
^max ^ 5 Q e y)_ p or eacn simulated test quark, its final location and momentum is sorted into this grid. To warrant 
a (statistically) sufficiently smooth behavior of the computed meson observables, a sample of ~10 8 test particles is 
needed. Finally, an interpolation algorithm has been devised for converting the discretized distribution back into a 
continuous function, to be plugged into the hadronization formula. The algorithm recovers the periodicity properties 
of the two angular variables, and converges to an arbitrary sampled function in the limit of a large number of grid 
points. A suitable grid dimension corresponding to the above variables amounts to, e.g., (13,96,14,12), where the 
large number of points in (f> p is dictated by the large sensitivity in the determination of the elliptic flow coefficient, 

V2{Pt)- 

We furthermore have to specify how to treat partons that escape the fireball prematurely, i.e., before the end of the 
QGP/mixed phase is reached. Clearly, these partons preferentially carry a high p t , and, upon exiting the fireball, could 
undergo a hadronization mechanism different from coalescence, such as fragmentation. Since our fireball is isotropic, 
the transition from the QGP to the vacuum is a sharp one and it would be unrealistic to coalesce the exiting quark 
with a thermal distribution at a temperature above T c , or at a time much later than the exit time (when the fireball 
has cooled down to T c ). We therefore decide to include only partons in our hadronization framework which remain 
inside the fireball throughout the entire QGP evolution. This leads to an underestimation of the high-momentum part 
of the hadronic spectra. A comprehensive calculation for quantitative comparison to experiment also at high p t would 
require the treatment of the exiting partons (hadronized with either fragmentation or coalescence). Similarly, the 



1 We recall [l2t | that the only requirement on the quark and meson masses is that the latter are above the two-quark threshold; within 
this restriction variations in the mass and (positive) Q values have little effect on the recombination process. In fact, as we will see 
below, CQNS scaling emerges approximately independent of quark mass. 
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FIG. 2: (Color online) Quark pt distributions resulting from relativistic Langevin simulations of an expanding elliptic QGP 
fireball at the end of a QGP (mixed) phase at a temperature of T = 180 MeV and average radial surface expansion velocity, 
(Po) — 0.48c. The numerically computed spectra (long-dashed lines) are compared to the initial spectra (solid lines) and to a 
blast-wave parameterization for particles with the same mass and the same fireball conditions (T = 180 MeV and the same flow 
field as used for the background medium in the Langevin simulation; short-dashed lines). Left and right panels correspond to 
charm (m c = 1.5 GeV) and strange (m s = 0.45 GeV) quarks, respectively. 



hadronic V2 we compute only reflects the partons within the QGP fireball at the end of its lifetime (which, however, do 
include non-thermal components from the Langevin simulation, in addition to the thermalized part of the spectrum). 

The generated spectrum, which originally represents a probability distribution, requires a suitable normalization. 
Since the empirical light and strange hadron spectra are consistent with chemical equilibrium close to the expected 
phase boundary, we assume this to apply at the quark level as well, at the critical temperature T c — 180 MeV of our 
fireball evolution. The fireball volume has been adjusted to match the total entropy of the fireball to the experimental 
hadronic final state multiplicities at T = 180 MeV at given collision centrality, e.g., Vfb — 1200 fm 3 for semicentral 
Au-Au collisions (note that one fireball covers approximately 1.8 units in rapidity). For charm quarks we augment the 
chemical equilibrium number by a fugacity factor, j c ~ 5 [331 ] . to match their number to the expected hard production 
in primordial nucleon-nucleon collisions (binary collision scaling; in Ref. [l2j 7c — 8 at T c = 170 MeV leads to the 
same number of cc pairs in the fireball). We note, however, that the overall normalization has little impact on our 
main considerations of space-momentum correlations and vi systematics. 

Let us finally specify the assumptions on the rapidity (y) distributions. As mentioned above, for both charm and 
strange quarks we employ a step function as 

S-SL* (*-«)■ 

where the parameter Ay depends on the quark mass and has been adjusted to recover approximately the full-width- 
half-maximum of a thermal y spectrum (amounting to Ay(s) = 1.3 and Ay(c) = 0.8 in connection with the quark 
masses quoted above). 

Fig. [^summarizes the results of the Langevin simulations for the p t probability distributions (integrated over spatial 
coordinates) for c and s quarks, compared to (normalized) blast-wave spectra for T = 180 MeV, using the same flow 
field as in the Langevin simulation. The average surface-expansion velocity is ((3q) ~ 0.48c. At low pt the spectra 
approach the equilibrium limit as to be expected since the background medium for the Langevin simulation is assumed 
to be fully equilibrated at all pt, as in a hydrodynamic calculation, and the stochastic process has been realized as to 
guarantee the correct equilibrium limit (including the adjustment of the (longitudinal) diffusion coefficient according 
to Einstein's fluctuation-dissipation relation [2lJ). As discussed above, the interaction strength has been chosen to 
recover the empirically observed maximum elliptic flow of v™ x ~ 7-8%. Note, however, that the assumption of a fully 
thermalized background medium implies rather large V2 values at high p t , while the phase-space density of thermal 
partons is rather small. A full treatment of this problem would require to solve a selfconsistcncy problem, where the 
parton spectra in the background fireball evolution also exhibit a saturation of the elliptic flow, ^(pt). This is beyond 
the scope of the present paper. 
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FIG. 3: (Color online) Meson pr spectra from quark-antiquark coalescence in central ^JIsnn = 200 GeV Au-Au collisions 
computed within the resonance recombination model for J/ip (left panel) and <f) (right panel); experimental data are from 
Refs. H and [13, HE HI, respectively. 



IV. MESON p T SPECTRA AND v 2 SYSTEM ATICS 



We now combine the ingredients of our resonance recombination model (RRM) by implementing the quark spectra 
computed from Langcvin simulations in the previous Section with the Boltzmann-based hadronization formalism 
of Sec. [Til evaluated in the stationary (equilibrium) limit according to Eq. ([7]). The key issue we address is how 
the properties of the input (non-observable) quark spectra reflect themselves in the (observable) meson spectra, in 
particular whether the space-momentum correlations generated in the Langcvin simulations can be consistent with 
the empirically observed CQNS, which, in turn, opens a window on the quark spectra at hadronization. 

It remains to specify the masses and widths of hadrons in the recombination process. In line with our restriction to 
mesons located above the quark-antiquark threshold (due to the limitation to 2 — > 1 processes) we consider s-s, c-c 
coalescence with resonance masses corresponding to the vacuum values for (j> (1.02 GeV) and J /if) (3.1 GeV) mesons; 
in connection with the quark masses as given above this implies similar Q values of 0.1-0.12 GeV. The (total) meson 
widths are chosen of comparable magnitude, i.e., = 0.05 GeV and Tj/^ = 0.1 GeV. As elaborated in Ref. [I4I, the 
numerical results, especially for the meson V2, are rather insensitive to variations in the meson width as long as Q is 
positive and substantially smaller (not smaller) than the resonance mass (width). 

In Fig. [3] we display our RRM results for pr spectra of J/ijj and cf> mesons, including available RHIC data. Overall, 
the spectra largely agree with those computed in our previous work (Fig. 1 in Ref. [l2l]). where the input spectra 
were solely based on a blast-wave parameterization. This is not surprising, since the quark spectra employed in the 
present work show a rather large degree of thcrmalization, up to momenta of pt = 1.5-2 GeV for strange and charm 
quarks, cf. Fig. [2j Consequently, the <p-m.eson spectra shown here are somewhat harder than in Ref. [l2| beyond 
Pt — 3 GeV due to the presence of the kinetic (hard) components in the quark spectra resulting from the Langevin 
evolution. The computed spectra for the J/tp are quite reminiscent to earlier blast- wave based results [321 . |34|. |35| . 
Note, however, that in Ref. (35[ the recombination (blast wave) contribution only amounts to about 50% of the total 
J /if) yield, significantly less than in the present paper (this is sensitive to the total open-charm cross section, which is 
not very well determined yet). 

The RRM results for the meson V2(pt) are summarized in Fig. [¥] (solid lines) and compared to the underlying 
quark scaled to meson variables in the conventional (empirical) way as i>2 Caled (pT') = 2^2,5 (pt/2). We find that for 
both the J /if) and <f> the agreement is rather impressive, within a few percent relative deviation. The wiggles at the 
quark level are, to a large extent, driven by the finite grid sampling due to a step width of 400 MeV (the statistical 
error inherent to the Langevin simulation is smaller than that, using 10 s test particles). The convolution of two 
quark distributions results in much smoother curves at the meson level (due to the fitting procedure of the quark 
input). We are thus able to approximately recover CQNS in a microscopic calculation with the full information on 
space-momentum correlations, characteristic for hydrodynamic expansion at low pt and a kinetic regime at higher p t . 

Another potential source for scaling violations are flavor (or mass) dependencies at the quark level (rather than 
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FIG. 5: (Color online) Elliptic flow coefficient, as a function of transverse momentum for c and s quarks (left panel), as well 
as J/ip and (f> mesons resulting from quark recombination (right panel), in semicentral ^/sjvjv = 200 GeV Au-Au collisions. 



in the coalescence process) . In Fig. [5] we compare the elliptic flow of different quarks (left panel) and mesons (right 
panel) with each other. The left panel confirms that the c-quark v-i deviates significantly (up to ~20% at low p t ) 
from the strange-quark V2- At high pt, the strange-quark V2 is a bit high, which, in principle, could be readjusted 
by a somewhat reduced strength of the resonance interactions in the Langevin simulations. At the meson level, the 
differences are similar. We have verified that comparable deviations persist when reducing the strange quark V2- 

Finally, we address the question of CQNS with respect to transverse kinetic energy, Kt = tiit — m, rather than 
transverse momentum, px- Such a scaling has recently been highlighted by the PHENIX [l9[ and STAR collabora- 
tions [2(| and seems to be very well satisfied by all available RHIC data, in centrality, collision energy and collision 
systems, after a geometric correction for the nuclear overlap. Kt scaling has also been found to result from cer- 
tain classes of hydrodynamic solutions , and therefore been argued to reflect a collectively expanding thermalized 
system of partons [41| . From the point of view of quark coalescence, the problem of reconciling quark distribution 
functions with space-momentum correlations (as implied by hydrodynamic expansion) with CQNS persists. In Fig. [5] 
we display the RRM results for V2, q and V2.m for the two different flavors as a function of K t ^r- We find that the quark 
input distributions from the QGP Langevin simulations indeed share a rather universal behavior up to K± ~ 3 GeV, 
encompassing both the quasi-equilibrium regime at low energies and the kinetic regime at intermediate energies char- 
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FIG. 6: (Color online) Elliptic flow coefficient, V2, as a function of transverse kinetic energy, Kt,T, for strange and charm quarks 
(left panel), as well as for <f> and J/tp mesons (right panel). 

acterized by a leveling off at K t > 1 GeV, cf. left panel of Fig. [6] We recall that the only adjusted input to this result 
is the common maximum value of the individual quark elliptic flow at about 7-8% (as suggested by the empirical 
CQNS deduced from experiment), controlled by the nonperturbative interaction strength in the stochastic Langevin 
force (again, a fine tuning for the s-quark would improve the agreement at higher K t ). The approximate universality 
at the quark level is nicely preserved at the meson level as a result of our Boltzmann-based recombination formalism, 
see right panel of Fig. [5] This is a quite remarkable result in view of the underlying space-momentum correlations in 
our approach, which has not been achieved before in this form. 



In the present paper we have extended our previously formulated quark coalescence formalism, utilizing resonance 
interactions within a Boltzmann equation, by implementing microscopic quark phase-space distributions generated 
via Langevin simulations of an expanding thermal QGP fireball for Au-Au collisions at RHIC. In this way we could 
combine the merits of our recombination approach (energy conservation and a proper equilibrium limit) with those 
of realistic quark distributions, which in particular encode the transition from a thermal regime at low pt to a 
kinetic one at intermediate p t . The latter feature is especially important as it produces the leveling-off of the elliptic 
flow, a key feature of observed hadron spectra and a crucial prerequisite to test any kind of quark scaling behavior. 
The (constituent) quark-mass and meson parameters were fixed at rather standard values, and we have constrained 
ourselves to mesons which arc reliably calculable in our 2—^1 recombination setup (i.e., for positive Q values). 
The only real adjustment concerned the interaction strength in the Langevin process, as to reproduce the empirical 
maximum value for the quark V2 (with the QGP fireball parameters tuned to empirical values of radial and elliptic 
flow, as in earlier applications to, e.g., heavy-quark observables). These interactions have been modeled via (meson) 
resonances in the QGP, which we identified with the states formed in the coalescence process at T c , leading to the 
notion of a "Resonance Recombination Model" (RRM). Since the Fokker-Planck approach as an expansion of the 
Boltzmann equation is strictly valid for sufficiently massive and/or high-momentum particles, we restricted ourselves 
to "heavy" flavors, i.e., charm and strange quarks. At low p t , the latter are at the borderline of applicability of a 
Fokker-Planck framework. 

Our main finding is that within this rather generic set-up, largely based on first principles augmented by a concrete 
realization of the strongly interacting QGP, the constituent quark scaling of the meson elliptic flow emerges rather 
naturally including space-momentum correlations characteristic for a collectively expanding source. The scaling 
holds for individual mesons, but appears to be rather universal in quark and meson flavor (mass), especially when 
applied in (transverse) kinetic energy rather than momentum, which is in line with recent experimental findings. By 
overcoming some of the limitations of previous (more schematic) coalescence models, and by achieving the first robust 
implementation of realistic (microscopically computed) phase-space distribution functions of quarks, our formalism 
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CONCLUSIONS 
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could provide a useful tool to better understand systematics of RHIC data, most notably the interplay of a thermal and 
kinetic regime in connection with phase-space properties of the partonic fireball as viewed through the hadronization 
process. 

Clearly, a formidable list of open issues persists, including the extension to light quarks, the role of gluons and 
of deeply bound hadronic states (possibly requiring additional formation processes), more realistic spectral functions 
of mesons and quarks and their interactions (both around T c and above) , a selfconsistent treatment of thermal and 
kinetic components (possibly requiring full parton transport), a systematic classification of viable parton phase-space 
distributions, hadronic reinteractions, etc. Progress has already been made on a number of these aspects, but a 
comprehensive approach remains a challenging task. 
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